Zero-field superfluid density in rf— wave superconductor evaluated from the results of 
muon-spin-rotation experiments in the mixed state 
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We report on measurements of the in-plane magnetic penetration Xab in the optimally doped 
cuprate superconductor (BiPb)2(SrLa)2Cu06+i (OP Bi2201) by means of muon-spin rotation 
(/iSR). We show that in unconventional d— wave superconductors (like OP Bi2201), nSK experi- 
ments conducted in various magnetic fields allow to evaluate the zero-field magnetic penetration 
depth Ao, which relates to the zero-field superfluid density in terms of ps oi X^^. 



PACS numbers: 74.72.Hs, 74.25. Jb, 76.75.-|-i 

Muon-spin-rotation (^SR) measurements in the mixed 
state of type-II superconductors provide valuable infor- 
mation on the superconducting properties. An impor- 
tant advantage of this method is that the muons probe 
the bulk of the material, and the results are not compli- 
cated by surface imperfections. The quantitative param- 
eters extracted from ^SR experiments depend, however, 
on the details of the model applied to reconstruct the in- 
ternal magnetic field distribution in the superconductor 
in the mixed state. So far. field distributions measured 
by means of /iSR were analyzed within the framework 
of analytical models based on London and Ginzburg- 
Landau (GL) theories, which can be applied, in general, 
to conventional superconductors with a single isotropic 
energy gappii^i^ The situation becomes much more com- 
plicated in the case of unconventional superconductors, 
like cuprates, MgB2, etc. It was found, in particular, that 
the effective magnetic field penetration depth Ae// ex- 
tracted from /LiSR measurements depends on the applied 
magnetic field ( see e.g. Refs. which is 

unexpected within the GL theory. Here Ae// refers to a 
quantity evaluated from ^SR experiments conducted in 
a superconductor in the mixed state, in contrast to Aq 
as obtained from Meissner state experiments {H ^ Hd, 
Hci denotes the lower critical field). In addition, it was 
observed that not only the absolute value, but also the 
shape of Kff{T) changes with field-^iii^ii^iii In this re- 
spect the question concerning the relation of Ag// to Aq, 
which is generally assumed to be proportional to the su- 
perfluid density (Aq ^ (x ps), becomes very important. 

In this paper we report on the results of a /LtSR study 
of the in-plane magnetic penetration depth Aab in op- 
timally doped (BiPb)2(SrLa)2Cu06+5. Ae//(T,iJ) was 
obtained from the measured temperature dependence of 
the /iSR linewidth by using numerical calculations of 
Brandt.-S. The temperature dependence of X^^ was fur- 
ther evaluated from Ae/ / (T, H) considering the nonlin- 
ear and the nonlocal response of a superconductor with 
nodes in the energy gap to the applied magnetic field. 
It was found that only at relatively low magnetic fields 
[B/Bc2{Q) ^ 10"'^, i3c2(0) is the zero-temperature value 



of the upper critical field] Ae// is a good measure of Aq. 
The high field data, however, need to be evaluated by 
taking into account both the nonlinear and the nonlocal 
corrections. 

Details on the sample preparation of optimally doped 
(BiPb)2(SrLa)2Cu06+5 (OP Bi2201) single crystals can 
be found elsewhere.— Field-cooled magnetization 
{Mpc) measurements of OP Bi2201 were performed with 
a SQUID magnetometer at ^qH — 1 mT, applied parallel 
to the c axis, for temperatures ranging from 5 K to 50 K. 
The transition temperature = 34.8 K was obtained as 
the intersect of the linearly extrapolated Mfc{T) curve 
in the vicinity of Tc with the M = line [see Fig. [1] (a)]. 
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FIG. 1: (Color online) (a) Field-cooled magnetization 
Mfc{T) of OP Bi2201. The field ^J.oH = 1 mT was ap- 
plied parallel to the crystallographic c axis, (b) The magnetic 
field distribution P{B) of OP Bi2201 taken at T = 1.6 K, 
fioH = 0.04 T. The lines represent the best fit within a two- 
Gaussian approach. 

The transverse-field /^SR experiments were carried out 
at the 7rM3 beam line at the Paul Scherrer Institute (Vil- 
ligen, Switzerland). Two OP Bi2201 single crystals with 
an approximate size of 4x2x0.1 mm'^ were mounted on 
a holder specially designed to perform /iSR experiments 
on thin single crystalline samples. The sample was field 
cooled from above Tc to 1.6 K in series of fields rang- 
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ing from 5 mT to 0.64 T. The magnetic field was ap- 
plied parallel to the c axis and transverse to the muon- 
spin polarization. The typical counting statistics were 
~ 15 — 18 million muon detections per data point. The 
experimental data were analyzed within the same scheme 
as described in Refs. iimiil. This is based on a two- 
component Gaussian fit of the /iSR time spectra which 
allows to describe the asymmetric local magnetic field 
distribution P{B) in the superconductor in the mixed 
state [see Fig. [1] (b)]. The magnetic field penetration 
depth A was derived from the second moment of P[B) 
as (T^ oc A~^ii The superconducting part of the square 
root of the second moment (usc oc ) was obtained 
by subtracting the normal state nuclear moment contri- 
bution {onm) from the measured cr, as = — cr^„ 
(see Ref. for details). Since the magnetic field was ap- 
plied along the crystallographic c axis, our experiments 
provide direct information on Aaf,. 

The temperature dependences of cTgc oc mea- 
sured after field-cooling the sample from far above in 
Aio-ff=0.04 T, 0.1 T, 0.2 T, 0.4 T, and 0.64 T are shown in 
Fig. [21 (a)- To ensure that a^dT) is determined primarily 
by the variance of the magnetic field due to the vortex 
lattice (VL) we plot in Fig. [2] (b) the corresponding val- 
ues of the skewness parameter as = (AB^)^^^ / {AB"^)^^^ 
[{AB") is the n— th central moment of P{B)]. as is a di- 
mensionless measure of the asymmetry of the lineshape, 
the variation of which reflects underlying changes in the 
vortex structureJ^ For an ideal triangular VL as ~ 1.2. 
It is very sensitive to structural changes of the VL which 
can occur as a function of temperature and/or magnetic 
Held^^ Fig.[2](b) implies that in OP Bi2201 as{T,H) 
is almost constant for 1.6 K< T < 26 K and smaller 
than the ideal value of 1.2, which is probably caused 
by distortions of the VL due to pinning effects. The 
sharp change of as at T ~ 30 K is similar to what 
was observed in Bi2212, where it was attributed to VL 
meltingj ^^i^^ Therefore, we conclude that for tempera- 
tures 1.6 K< T < 30 K the T variation of (Tsc in OP 
Bi2201 studied in the present work reflects the intrin- 
sic behavior of the in-plane magnetic penetration depth 

\ab{T). 

From the measured asc{T,H) we reconstructed 
A~j J (T, H) by using the procedure described in Ref. [l3. 
A correction between dsc and A~y^y was considered ac- 
cording to: 



,{h)[^s-^]= A{b)\ 



(1) 



which accounts for decreasing of the field variance within 
the VL with increasing magnetic fieldjii^ The correction 
factor A{b) depends only on the reduced field b = B/Bc2 
{Bc2 is the upper critical field). For a superconductor 
with a Ginzburg-Landau parameter k = A/f > 5 mea- 
sured in fields ranging from 0.25/ K^^^ ^ b < 1, A{b) can 
be obtained analytically as A{b) = 4.83 • 10'*(1 - &)[1 + 
1.21(1 - Vb)^] fis-^nm^ (see Ref.i). 

The reconstructed \~^^{T^H = const) curves are 
shown in Fig. [3l The corresponding A(6) dependences 
are displayed in the inset. The calculations were made 
for i?c2(0) = 50 Tfi^ The temperature dependence 




FIG. 2: (Color online) (a) Temperature dependence of dsc oc 
A~^^ of OP Bi2201 measured at [loH = 0.04 T, 0.1 T, 0.2 T, 
0.4 T, and 0.64 T. (b) Dependence of the skewness parameter 
Qfs on temperature. 
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FIG. 3: (Color online) Temperature dependence of the ef- 
fective magnetic penetration depth reconstructed from 
asc{T) measured at [iqK = 0.04 T, 0.1 T, 0.2 T, 0.4 T, and 
0.64 T (see Fig. [Sjl. The inset shows the temperature depen- 
dence of the correction factor A{b) — asc ■ ^eff- 



of Bc2 was assumed to follow the Werthamer-Helfand- 



Hohenberg (WHH) prediction^ Below T - 20 K, A 
is linear in T, as is expected for superconductor with 
nodes in the energy gap. Fig. [3] also implies that in 
the whole temperature region (from T ~ 1.6 K up to 
Tc), \~^j:{T, H) decreases with increasing field. This 
contrasts the results obtained by using a similar pro- 
cedure for the ternary boride Li2Pd3B and electron- 
doped Sro.9Lao.iCu02i ^^i^° For these two compounds the 
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(T, H) curves were found to collapse onto a single 
curve. Since Li2Pcl3B and Sro.9Lao.iCu02 are supposed 
to be fully gaped )^''i^°'^^i^^'^'^i^^ we may conclude that the 
field dependence of \~J^{T), shown in Fig. [31 is caused 
by the presence of nodes in the superconducting energy 
gap of OP Bi2201,25 
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FIG. 4: (Color online) Dependence of the proportionality fac- 
tor C(fe) = Ae/ / (T, b) / Ao (T) on the reduced field 6 = B /B^2 of 
OP Bi2201. The solid Une is the fit by means of Eq. © which 
takes into account the nonlinear correction to Aq. The inset 
shows the magnetic field dependence of \eff aX T = 1.6 K, 
10 K, 17 K, and 23 K. The reduced fields for two points se- 
lected in the oval are almost the same. 

As shown in Refs. [2^ and [l^, the magnetic field de- 
pendence of Ae// arises from the nonlocal and the non- 
linear response of a superconductor with nodes in the 
energy gap to the applied magnetic field. The nonlinear 
correction to Aab appears due to the magnetic field in- 
duced quasiparticle excitation over the gap nodes.— Ac- 
cording to Volovik,— the density of the dclocalizcd states 
increases proportionally to ^/h. The nonlocal correction 
to \ab appears from the response of electrons with mo- 
menta on the Fermi surface close to the gap nodes. This 
is because the coherence length ^, being inversely propor- 
tional to the gap, becomes very large close to the nodes 
and, formally, diverges at the nodal points. Thus there 
exist areas on the Fermi surface where A/^ < 1, and the 
response of a superconductor to an applied magnetic field 
becomes highly nonlocali^ 

In order to reconstruct the temperature dependence 
of the superfluid density in zero magnetic field we used 
the general assumption that the proportionality factor 
relating Ae/ / to Aq is a function of the reduced magnetic 
field b only, so that: 



E 4 




1.0 



s o.i 



0.6 



0.4 



0.2 



^ 0.0 



(b) 



X aJ0.04T) 

* \„"(0.04T) 

□ 0^{O.64T) 

. i,/(0.64T) 

10 



20 

r(K) 



30 



FIG. 5: (Color online) (a) Temperature dependence of Aq ^ 
of OP Bi2201 reconstructed from \~jj{T,H) measured at 
/ioi^ = 0.04 T, 0.1 T, 0.2 T, 0.4 T, and 0.64 T (see 
Fig. El . (b) Temperature dependence of the normalized X^^ , 
A .V0.04 T), cr,,(0.04 T), Aj' (0.64 T), and a,, (0.64 T). 



selected temperatures. It is seen, e.g., that the reduced 
field b for Ae// measured at T = 10 K, iiqH = 0.4 T 
(lower point in the oval selection) is almost the same as 
the one for the point at T = 23 K, iiqH ~ 0.2 T (upper 
point). This implies that the coefficients C(6) for these 
two points are nearly equal and that the difference in 
the absolute values of Ae//(10 K, 0.4 T) and Ae//(23 K, 
0.2 T) is due to different values of Aq. The reconstruc- 
tion procedure was performed in the following way. First, 
from \-ff{T,H) plotted in Fig. [3 Ae//(5) was recon- 
structed for various constant temperatures (see inset in 
Fig. [4]). Second, the resulting values of Ae//(T = const, b) 
were scaled in order to have them collapsing on a single 
curve (see Fig. According to Eq. ^ this curve cor- 
responds to C{b) — Xef f{T,b)/ Xo{T), while the scaling 
factor, in turn, corresponds to Ao(r). The solid line rep- 
resents the result of the fit by means of the relation: 



Ae//(6,T) = C(6)Ao(T). 



(2) 



Ae//W/Ao = C{b) = (1 - KVb)-^/\ 



(3) 



This statement is correct, at least, in case of nonlinear 
corrections which scale with ^^, 28, 29^30 

The fact that the field and the temperature depen- 
dences of Ae// are described by separate terms [see 
Eq. ([2])] allows to reconstruct Ao(T'). In order to demon- 
strate this, we refer to the inset in Fig. 0] which shows 
the dependence of Ae// on the magnetic field for some 



which takes into account the nonlinear correction to Aq 
for a superconductor with d— wave energy gap.— >^ Here 
the parameter K depends on the strength of the non- 
linear effect. It is obvious that the "nonlinear" curve 
describes the experimental C{b) = Xeff{b)/XQ depen- 
dence reasonably well. In particular, it reproduces the 
curvature at 5 < 0.01 and the linear increase of C{b) for 



4 



0.01 ^ & ^ 0.05. We believe, however, that the whole 
Ae//(&)/Ao curve must be a combination of both nonlin- 
ear and nonlocal correction effects, similar to the results 
of Ref . m. 

The temperature dependence of Aq ^ and the compar- 
ison of Xq^{T) with (Tsc{T) and \~J^-{T) measured at 
Moi? = 0.04 T and 0.64 T arc presented in Fig. [H Both 
asc{T) and \~^j{T) measured at ^qH = 0.04 T almost 

coincide with each other as well as with AQ'^(r). This 
implies that the two sets of corrections, namely, the first, 
accounting for decrease of the second moment of the /xSR 
line with increasing field [Eq. ^ and the inset in Fig. [5] 
and the second, arising due to the nonlocal and the non- 
linear response of a superconductor with nodes in the 
gap to the apphed magnetic field [Eq. ^ and Fig.|3], are 
not really important at this relatively low field. Conse- 
quently, the second moment of /iSR line ct^c measured at 
/io-ff = 0.04 T is still a good measure of the zero-field su- 
perfluid density ps oc A,^^. On the other hand, asc{T) and 
\~ff{T) bX fioH ~ 0.64 T differ substantially from each 

other and from the resulting AQ^(r). This implies that 
for high fields both above mentioned corrections have to 



be taken into account. 

To conclude, muon-spin rotation measurements were 
performed on the optimally doped cuprate supercon- 
ductor (BiPb)2(SrLa)2Cu06+5. It was demonstrated 
that in unconventional d— wave superconductors (like OP 
Bi2201) /xSR experiments taken in various magnetic fields 
allow a reliable evaluation of the zero-field supcrfluid 
density ps oc Aq"^. The deviation of the effective pen- 
etration depth Ae// from Aq observed for higher fields 
was explained by the nonlinear and the nonlocal re- 
sponse of the superconductor with nodes in the energy 
gap to the applied magnetic field. The dependence of 
Ae///Ao on the reduced magnetic field h = B /Bc2 follows 
a (1 — K^/b)~^/^ behavior, accounting for the nonlinear 
correction to Aq for a wave superconductor. 
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